Relativistic plasma control for single attosecond pulse generation 
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To describe the high harmonic generation at plasma surfaces in the relativistic regime, we introduce 
the concept of apparent reflection point (ARP). It appears to an external observer that the radiation 
electric field is zero at the ARP. The relativistic dynamics of the ARP completely defines the 
generation of high harmonics and attosecond pulses. The ARP velocity is a smooth function of 
time. The corresponding 7-factor, however, has sharp spikes at the times when the tangential 
vector potential vanishes and the surface velocity becomes close to the speed of light. We show 
that managing the laser polarization, one can efficiently control the ARP dynamics, e.g., to gate a 
single (sub-) attosecond pulse out of the short pulse train generated by a multi-cycle driver. This 
relativistic control is demonstrated numerically by particle-in-cell simulations. 

PACS numbers: 52.27.Ny,42.65.Ky 



High harmonics generated at plasma surfaces in the 
relativistic regime []]■ El are a promising new source of 
short wavelength radiation and attosecond pulses, nec- 
essary for the study of ultra- fast processes in atoms, 
molecules and solids at intensities significantly higher 
than those obtained from strong field laser-atom inter- 
actions nm. 

The plasma harmonics are generated in the relativistic 
regime, when the laser pulse intensity / 3> 10 18 W/cm 2 . 
Laser pulses in this intensity range usually are several 
cycles long. As a consequence, the reflected radiation 
contains a comb of attosecond pulses. Yet, applications 
like molecular imaging 0, [fj or quantum control usu- 
ally require a single short pulse to prevent undesirable 
effects, e.g., Coulomb explosion. The single attosecond 
pulse can be selected either using a phase-stabilized sin- 
gle cycle laser pj or controlling the atomic response by 
time-dependent laser polarization 0, || . 

It was shown recently by particle-in-cell (PIC) simu- 
lations that in the A 3 — regime, when a single-cycle laser 
pulse is focused down to a spot of a one wavelength size, 
a single attosecond pul se running in a particular direc- 
tion can be isolated [jjj]. However, in the case of a more 
realistic several-cycle laser pulse, such casual separation 
can be difficult. 

In the present work we show for the first time that 
the managed time-dependent polarization (from elliptical 
to linear and back) of the incident laser pulse allows to 
gate a single (sub-)attosecond pulse from relativistically 
driven plasma surface in a well controlled way. Physi- 
cally, the laser pulse polarization controls the relativistic 
7— factor of the apparent reflection point (ARP) as seen 
by the observer. Although the driving laser pulse can be 
long and intense, the ARP velocity becomes highly rela- 
tivistic only at the specific time, when the vector poten- 
tial component tangential to the plasma surface vanishes. 
Exactly at this moment, the single attosecond pulse is 
emitted. The time-dependent polarization corresponds 
to two perpendicularly polarized laser pulses with slightly 
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FIG. 1: Geometry for relativistic plasma control of attosecond 
plasma surface dynamics. The high intensity driving pulse 
is polarization-managed through a low intensity controlling 
pulse. After proper filtering of the reflected radiation, a single 
attosecond pulse can be isolated. The interaction process can 
be described in terms of the apparent reflection point (ARP). 
It appears to an external observer that the radiation electric 
field and the Poynting vector S are zero at the ARP. 



different frequencies and a well chosen phase shift, see 
Fig. Experimentally, the time-dependent laser pulse 
cllipticity can be achieved by the femtosecond polariza- 
tion pulse shaping techniques 0. 

We consider the interaction of an ultra-intense short 
laser pulse with a slab of overdense plasma. We suppose 
that the plasma ions are immobile during the short inter- 
action time and study the electron fluid oscillations only. 
It was shown in |3( that to describe the high harmonic 
generation analytically, the boundary condition E r = 
must be used. Here E r is the electric field component 
tangential to the plasma surface. The physical meaning 
of this boundary condition is quite clear. Let us con- 
sider the Poynting vector S = cE x B/47T that gives 
the electromagnetic energy flux. In the vacuum region 
in front of the plasma, it consists of the two parts: the 
incident flux S 1 = cE 1 x H 1 /4tt and the reflected one 
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S r = cE r x B r /47r. Because the plasma slab is overdense, 
S = behind the plasma. If we also neglect the small ab- 
sorption, then the observer sees the apparent reflection 
at the point x AKP (t), where the normal component of 
the Poynting vector S n = cE T x B r /47r = 0, implied by 
E T (x ARP ) = 0. Thus the boundary condition ensures the 
energy conservation. The detailed microscopic derivation 
of the boundary condition E T (x ARP ) = and its relation 
to the electron dynamics can be found in ^1 ■ 

The incident laser field in vacuum runs in the positive 
x— direction, W{x,t) = W{x — ct), while the reflected 
field is translated backwards: W(x,t) — ~E r (x + ct). The 
tangential components of these fields interfere destruc- 
tively at the ARP position x AKP (t), so that the implicit 
equation for the apparent reflection point x AB _ P (t) is: 

E^(z ARP - ct) + K(xaup + ct) = 0. (1) 

We stress that Eq. Q contains the electromagnetic fields 
in vacuum. That is why the reflection point x ARr is ap- 
parent. The real interaction within the plasma skin layer 
can be very complex. Yet, an external observer, who 
has information about the radiation in vacuum only, sees 
that E r = at a; ARP . Strictly speaking, one cannot as- 
cribe a: ARP to any particular density level, e.g., the "criti- 
cal density" . One can only state that the ARP is located 
within the skin layer at the electron fluid surface. 

The ARP dynamics completely defines the high har- 
monic generation. In the ultra-relativistic regime, when 
the dimensionless vector potential ao = eA^/mc 2 of the 
laser is large, al_^> 1, we can apply ultra-relativistic sim- 
ilarity theory [lj to characterize this motion. The ba- 
sic statement of this theory is that when we change the 
plasma density N e and the laser amplitude ao simultane- 
ously keeping the similarity parameter S — N e /aoN c con- 
stant the laser-plasma dynamics remains similar. Here 
N c = Lo^m/Ane 2 is the critical plasma density for the 
laser pulse with the frequency luq. This means that for 
different interactions with the same similarity parameter 
S'=const, the plasma electrons move along the same tra- 
jectories, while their momenta p scale with the laser am- 
plitude: p rx ao- Consequently, the electron momentum 
components tangential and normal to the plasma sur- 
face scale simultaneously with ao: p T oc ao and p n cx ao- 
Therefore, in general, the total electron momenta are not 
perpendicular to the surface. Moreover, the characteris- 
tic angle between their direction and the surface normal 
does not depend on ao provided that S is fixed. 

This result is crucial for the plasma surface dynamics. 
Since we consider an ultra-relativistic laser pulse, a\ ^> 1, 
the electrons in the skin layer move with ultra-relativistic 
velocities almost all the time: 
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Yet the relativistic 7-factor of the plasma surface 7 S (t) 
and its velocity v s (t) behave in a quite different way. Let 



us consider electrons at the very boundary of the plasma. 
The similarity theory states that the electron momenta 
can be represented as p n (t) = a P n {S, wt) and p r (t) = 
aoP r (S, cut), where P n and P T are universal functions, 
which do not depend on ao explicitly, but rather on the 
similarity parameter S and on the pulse shape. Thus, 
the plasma surface velocity (3 S (t) — v s (t) /c and 7 S (t) are 
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It follows from Q-0 that the relativistic 7-factor of 
the plasma boundary is of the order of unity for almost 
all times, except for those times t g , when the tangential 
momentum component vanishes p T (S,t g ) = 0. Exactly 
at these times, there are spikes of the 7-factor: 



'Pn + mic 



2 r 2 



oc ao. 



(5) 



On the contrary, the plasma surface velocity v s as given 
by Eq. is a smooth function. It approaches ±c when 
the electron momentum parallel to the surface vanishes, 
i.e., at the same times t g corresponding to the 7-spikes. 
The high harmonics are generated at the spikes, when the 
surface velocity is negative and close to — c. Rigorous mi- 
croscopic analysis confirms that high harmonics are 
generated by bunches of fast electrons moving towards 
the laser pulse. 

It is possible to estimate the width of a 7-spike. Since 
the surface velocity is a smooth function, we can expand 
it in Taylor series around the maximum as v s (t) w v max — 
au 2 (t — t max ) 2 , where the parameter a depends only on 
S. The width At = \t — t max \ of the 7 S spikes is then 



Ai « V 1 - v max/ c2 /i UJ oVa) = 1/(^0 Va^fmax) ■ (6) 

It follows from @ that the 7 S spikes get higher and nar- 
rower when we increase ao, keeping S = const. Thus, 
the spikes behave as quasi-singularities. These ultra- 
relativistic spikes are the inherent cause for the high har- 
monic generation. The (sub-)attosecond pulses |3( are 
emitted exactly at the times of the spikes. 

We study the motion of the plasma boundary and 
the specific behavior of w ARP and 7 ARP numerically us- 
ing the ID PIC code VLPL The plasma slab 
is initially positioned between xl = 1.5A and xr = 
3.9A, where A = 27r/wo is the laser wavelength. The 
laser pulse has the gaussian envelope: a(x, t = 0) = 
ao exp(— x 2 /a 2 ) cos(27rx/A) with a — 2A. 
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FIG. 2: ID PIC simulation results for the parameters ao = 20 
and A^ e = 90^. a) Oscillatory motion of the point a; A Rp(t) 
where E T (x(t)) = 0. b) Velocity w A Rp(i) = dx a-rp (t) / dt; only 
the negative velocities are shown. Notice that the ARP ve- 
locity is a smooth function, c) The corresponding 7— factor 
7arp(£) = — «ARp(i) 2 /c 2 contains sharp spikes, which 



coincide with the velocity extrema. 



At every time step, the incident and the reflected fields 
are recorded at x — 0. Being solutions of the wave equa- 
tion in vacuum, these fields can be easily chased to arbi- 
trary x and t. To find the ARP position x ARP , we solve 
numerically the equation l|Tjl. The trajectory of x ARP (i) 
for the simulation parameters ao = 20 and N e /N c = 90 
(S = 4.5) is presented in Fig. [21 a). One can clearly see 
the oscillatory motion of the point x AnF (t). The equilib- 
rium position is displaced from the initial plasma bound- 
ary position xl due to the mean laser light pressure. 

Since only the ARP motion towards the laser pulse is 
of importance for the high harmonic generation, we cut 
out the positive ARP velocities v AUP (t) = dx ARP (t)/dt 
and calculate only the negative ones, Fig. [21 b). The 
corresponding 7- factor 7 ARP (t) — 1/ \/l — v AUP (t) 2 /c 2 



presented in Fig. [21 c) . Notice that the ARP velocity 
is a smooth function. At the same time, the 7-factor 
7 ARP (i) contains sharp spikes, which coincide with the 
velocity extrema. These spikes of the surface 7-factor 
are responsible for the high harmonic generation. 
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FIG. 3: a) The high harmonics spectrum; b) the train of 
attosecond pulses is obtained after the proper filtering. 



The numerically obtained spectrum of the high- 
harmonics is shown in Fig. [3| a). Filtering out the 
lower harmonics and keeping only the harmonics with 
lo > 25loq, we obtain a train of short pulses in the re- 
flected radiation, Fig. [31 b). 

For various applications, such as molecule imaging, it is 
of great importance to have a single short pulse, instead 
of a train of short pulses. We have shown above that 
the attosecond pulses are emitted when the tangential 
components of the surface electron momentum vanish. 
This property can be used to control the high harmonic 
generation and to gate a particular attosecond pulse out 
of the train. 

In the ID geometry, the transverse generalized momen- 
tum is conserved: p T = eA T /c, where p T and A T are the 
tangential components of the electron momentum p and 
the vector potential A. Consequently, the attosecond 
pulses are emitted when the vector potential is zero. If 
the vector potential vanishes at several moments, there 
are several 7-spikes and correspondingly, after proper fil- 
tering, several short pulses are observed in the reflected 
radiation, see Fig. [31 b). 

To select a single attosecond pulse, we must ensure that 
the vector potential A r turns zero exactly once. Since A r 
has two components, how often it vanishes depends on 
its polarization. For linear polarization it vanishes twice 
per laser period, while for circular or elliptic polarization 
it never equals zero. A laser pulse with a time-dependent 
polarization can be prepared in such a way that its vector 
potential turns zero just once. The time-dependent po- 
larization corresponds to two perpendicularly polarized 
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FIG. 4: Generation of a single attosecond pulse using the 
relativistic plasma control. The driver signal has ao = 20 
and frequency u>q = 1. The controlling signal has a c = 6 and 
frequency uj c = 1.25. The phase difference is A<f) — n/8. a) 
the 7-spikes of the oscillating ARP, b) the single attosecond 
pulse selected via the relativistic plasma control. 



laser pulses with slightly different frequencies, see Fig.^ 
Our PIC simulations suggest that a signal with a few per- 
cent of the driver intensity is sufficient to control the high 
harmonic generation, if the phase difference between the 
two laser pulses is chosen carefully. 

To demonstrate the relativistic plasma control, we per- 



form a PIC simulation where we add the z-polarized 
controlling pulse with amplitude a c = 6 and frequency 
oj c = 1.25 and retain the same driver pulse with ampli- 
tude ao = 20, frequency luq = 1 and y-polarization. We 
keep the same plasma density N e /N c — 90. The op- 
timal phase difference between the two lasers is found 
empirically to be Acf> — ir/8. The simulation results are 
presented in Fig. 0] a) and b) . Comparing the surface 
7— factor dynamics in the regime of linear polarization, 
Fig. |21 c) and in the controlled regime, Fig. 01 a), we see 
that the central 7— spike is slightly larger while the both 
side spikes are significantly damped. This effect becomes 
much more pronounced when we compare the filtered ra- 
diation plots, Fig. |21 b) and Fig. 01 b). The control signal 
allows us to select the single attosecond pulse correspond- 
ing to the highest 7— spike in the surface motion. 

Varying the control parameters ao/a c , ujq/lo c and A</> 
we are able to select different attosecond pulses one-by- 
one or in groups out of the original pulse train. 

To recapitulate, we studied analytically and numeri- 
cally the dynamics of the apparent reflection point at 
the overdense plasma surface. We have shown that the 
velocity of this point is a smooth function of time. How- 
ever, the corresponding 7-factor has quasi-singularities 
or spikes when the surface velocity approaches the speed 
of light. These ultra-relativistic spikes are responsible for 
the high harmonic generation in the form of an attosec- 
ond pulse train. We show that the attosecond pulse emis- 
sion can be efficiently controlled by managing the laser 
polarization. This is done by adding a low intensity con- 
trol pulse with perpendicular polarization and frequency 
slightly different from that of the driving pulse. This rel- 
ativistic plasma control allows to gate a single attosecond 
pulse or a prescribed group of attosecond pulses. 

This work has been supported in parts by DFG Tran- 
sregio 18 and by DFG Graduierten Kolleg 1203. 
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